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Electroweak baryogenesis (EWBG) in extensions of the Standard Model will be 
tested quantitatively in upcoming nuclear and particle physics experiments, but only 
to the extent that theoretical computations are robust. Currently there exist orders- 
of-magnitude discrepancies between treatments of charge transport dynamics during 
EWBG performed by different groups, each relying on different sets of approxima- 
tions. In this work, we introduce a consistent power counting scheme (in ratios of 
length scales) for treating systematically the dynamics of EWBG: CP-asymmetric 
flavor oscillations, collisions, and diffusion. Within the context of a simplified model 
of EWBG, we derive the relevant Boltzmann equations using non-equilibrium field 
theory, and solve them exactly without ansatz for the functional form of the density 
matrices. We demonstrate the existence of a resonant enhancement in charge pro- 
duction when the flavor oscillation length is comparable to the wall thickness. We 
compare our results with the existing treatment of EWBG by Konstandin, Prokopec, 
Schmidt, and Seco (KPSS) who previously identified the importance of flavor oscil- 
lations in EWBG. We conclude: (i) the power counting of KPSS breaks down in 
the resonant regime, and (ii) this leads to substantial underestimation of the charge 
generated in the unbroken phase, and potentially of the final baryon asymmetry. 
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I. INTRODUCTION 

Electroweak baryogenesis (EWBG) is an attractive mechanism to explain the origin of 
the baryon asymmetry of the Universe (BAU). Because the relevant physics is accessible at 
the electroweak scale, nuclear and particle physics experiments can directly test two of the 
three Sakharov conditions [1] . The departure from equilibrium during the electroweak phase 
transition (EWPT) can be addressed in collider studies (e.g., [2]), while CP violation can 
be probed through low-energy observables such as electric dipole moments (EDM) [3]. 

In EWBG, electroweak symmetry breaking proceeds through a first-order EWPT, where 
bubbles of broken SU(2)ixU(l)y symmetry nucleate and expand in a background of unbro- 
ken symmetry. CP-violating interactions within the bubble wall produce a CP-asymmetric 
charge density of left-handed fermions. This charge, diffusing ahead of the bubble wall into 
the unbroken phase, is converted into a baryon asymmetry through non-perturbative elec- 
troweak sphaleron processes [3]. Lastly, the baryon asymmetry is captured by the advancing 
bubble wall and freezes out. The Sakharov conditions are satisfied provided (i) the EWPT 
is "strongly" first-order, otherwise electroweak sphalerons are too active in the broken phase 
and the baryon density is washed out, and (ii) there exists sufficient CP violation to gener- 
ate the observed BAU. Although neither condition is met in the Standard Model (SM) [5], 
EWBG may be viable in the Minimal Supersymmetric SM (MSSM) [31 16] and other scenarios 
beyond the SM 0. 

Ultimately, the final baryon asymmetry in EWBG is roughly proportional to the total CP- 
asymmetric charge that is generated and transported, by diffusion, into the unbroken phase 
where weak sphalerons are active. However, within a given model, this charge transport 
computation is a complicated problem: one must solve a network of Boltzmann equations 
governing the generation, diffusion, and equilibration of charges in the vicinity of the moving 
bubble wall. To date, there exists no treatment of EWBG dynamics that includes all of these 
aspects in a consistent and theoretically robust framework. And yet, such a treatment is 
crucial for making a quantitative connection between EWBG and experiment. Here, a key 
insight was provided by Refs. [HI HD] in recognizing the importance of CP-violating flavor 
oscillations. These oscillations arise through spacetime-dependent flavor mixing induced by 
the varying background Higgs field. 

In the thick bubble wall regime (L w ^> T _1 ), flavor oscillations are formally the leading 
source of CP violation in a gradient expansion in powers of (L w T)~ l , arising at linear 
order (e.g., in the MSSM, L w ~ 20 /T |8J). At order (Z^X 1 ) ^ . one finds an additional 
CP-violating source from the spin-dependent "semi-classical force" [TTHT5] . Aside from 
Refs. O [10] , all previous EWBG computations have swept flavor oscillations under the rug: 
quantum coherence between states was neglected, and this CP-violating source was treated 
perturbatively in terms of scattering off the external background Higgs field or its gradient. 
It is unknown whether or not such prescriptions are justified. Furthermore, Ref. [H] found 
that CP-violating charge densities generated by flavor mixing are localized to the bubble 
wall and do not diffuse into the unbroken phase, in contrast to other treatments. On the 
other hand, if significant charge diffusion occurs it would lead to a substantial enhancement 
of baryon number generation by transporting charge outside the bubble where sphalerons 
are active [25T - T27] . Therefore, this issue is highly relevant for experimental tests of EWBG. 
Indeed, according to Ref. [9], EWBG in the MSSM is ruled out due to null electron and 
neutron electric dipole moment (EDM) searches, while it is still viable (to varying degrees) 
according to other EWBG computations [31 E]- 
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It is clear that quantitative EWBG computations require a consistent analysis of CP- 
violating flavor oscillations, collisions in the plasma, and diffusion, which is still missing in 
the literature. In Ref. [TB], we performed a first step towards such a complete treatment: we 
studied a model of two scalars fields $l,r that mix through a purely ttme-dependent CP- 
violating mass matrix and have flavor-sensitive interactions with a thermal bath of scalar 
bosons (A) in equilibrium. We derived using non-equilibrium field theory [17] quantum 
Boltzmann equations for the two-flavor density matrices (Wigner functions), relying on a 
systematic expansion scheme in ratios of time scales. Our work demonstrated — for the first 
time in a consistent framework — the resonant baryogenesis regime discussed in Refs. [18]- 
21] and placed it on a more rigorous theoretical footing. Importantly, we showed that the 
resonance occurs when the flavor oscillation period is comparable to the variation scale of 
the background Higgs field. 

In this work, we generalize the analysis of Ref. P5] to account for a spacetime-dependent 
background field geometry, as in a moving bubble wall. Non- homogeneity of the bubble 
geometry is a key ingredient of EWBG, essential for generation of a CP-violating charge 
that undergoes subsequent transport by diffusion. We work within the context of the same 
toy model of Ref. [16]. Our work is organized as follows: 



In Sec. [ITJ we present our toy model for two mixing scalars (&l,r) with a spacetime- 
dependent mass matrix induced by the bubble wall. We discuss how this model maps 
onto more "realistic" baryogenesis models. 

We derive the Boltzmann equations for the occupation numbers and quantum coher- 



ence of &l,r states in Sec. Ill We present both a heuristic, intuitive derivation, and 
a more rigorous one using non-equilibrium field theory. We also discuss the necessary 
conditions for CP violation. 



• Using numerical methods, we solve the Boltzmann equations exactly in Sec. IV, with- 
out ansatz for the form of the density matrices. Our results clearly demonstrate the 
existence of charge diffusion, the role of flavor oscillations in generating CP asymme- 
tries, and how charge generation is enhanced in the resonant regime \m\— m2\ < 10L" 1 , 
where are the mass eigenvalues of the two-scalar system. Here, we also provide a 
useful analogy with spin precession in a varying magnetic field. 

• In Sec. |VJ we highlight the differences between our approach and that of Refs. [HJ ITU] . 
Our major disagreement stems from Refs. [HI EE] power counting the solutions to the 
Boltzmann equations in powers of (L^T) -1 , while we do not. We demonstrate that 
this power-counting argument breaks down in the resonant regime, effectively negates 
the possibility of diffusion, and substantially underestimates the amount of charge 
generated during the EWPT. 

The virtue of working within our toy model is that it can be solved both exactly (with 
numerical techniques), as well as within various approximation schemes. This offers the 
possibility to study some of the key assumptions used in current approaches and to quantify 
the attendant uncertainties. In this work we have focused on the comparison with what has 
so far been considered the state-of-the-art calculation in Ref. [9j. In forthcoming work, we 
will study in detail the diffusion approximation, invoked in essentially all EWBG calculations 
to make the problem tractable. In future work we will also extend the application of our 
methods to fermions, which are an essential ingredient in EWBG since sphalerons couple 
only to fermions. 
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II. BARYOGENESIS TOY MODEL 



Two-flavor dynamics of scalars are highly relevant in extensions of the MSSM, where top 
squarks (t L , t R ) may account for the BAU [33J. These models are necessarily complicated by 
their large number of degrees of freedom ~ 200) and many different types of interactions. 
Here, we consider a much simplified model: a two-flavor scalar system, with fields $ = 
<3>#), described by the Lagrangian 

& = fy$t _ $t M 2 $ + ^f. nt . (!) 

The key ingredients of our model are: 

• The mass matrix M 2 (x) is spacetime-dependent, assumed to a function of the varying 
background Higgs field(s) associated with the expanding bubble. The variation of M 2 
across the phase boundary generates $>l,r charge through CP-violating coherent flavor 
oscillations. 



• J2?i nt describes the interactions of $ with the remaining degrees of freedom in the 
plasma. These collisions govern the damping of flavor oscillations and the effective 
diffusion of these charges in the plasma. 

We model the plasma during the time of the EWPT as a thermal bath of real scalar bosons 
A, assumed to be in equilibrium at temperature T, coupled to $ via 

^int=-^ 2 $t y$; „=f«f M. (2) 



2 " ' y V VR. 

We take the matrix of coupling constants y to be diagonal; this defines the basis of flavor 
eigenstate fields <3?l,r ("flavor basis"). 

The mass matrix can be parametrized as 

M 2 (t \ _ ( m l( x ) v ( x ) e ' ia(x) \ (o\ 
M [X) ~ \ v(x) e ia ^ rn 2 R (x) J ' [6) 

It is convenient to transform Eq. ([I]) into the basis of local mass eigenstates ("mass basis"). 
We diagonalize the mass matrix with the spacetime-dependent transformation matrix U(x), 
such that 

2/ \ f rn 2 (x) \ r7+ ., 2 rr tti \ f cos9(x) — sm6(x) e~ ia ^\ , 
m 2 (x )=[ ^ ; 2{ \ )=WM 2 U, U{x)=[ . a( ^ L x) K J, v 4 

v ' \ m^x) J y ' \sm9(x) e ia{x> cos9(x) J y ' 

with 

I/O 9\,1. / O o \ fl n n \ r, . '. o /„ ^\ 2 V 

'/. " "<lt 



m l,2 = ~( m i + m R) ± ^sign(m| - m 2 R )J{m 2 L - m 2 R ) 2 + 4v 2 , tan(20) = 2 _ 2 ■ (5) 



This diagonalization defines the mass basis fields = 02) = W$>. The Lagrangian, in 
the mass basis, is 

S? = <9 M 0t - - 0t E m 9 ^ + 9 S m _ S m E ^ + ^ ; (6) 
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where 

E» W S U\ X ) rTU {x) = m + (j ^ * *■« . (7) 
In this basis, the interaction becomes 

^int = - ^A 2 0t Y , y (x) = U\x) y U(x) . (8) 

During the EWPT, the background Higgs field has a bubble geometry. We assume a 
spherical bubble expanding in the r direction, with wall thickness L w ^> T _1 and velocity 
v w < 1. Typically, in the MSSM, one finds L w ~ 20/T [8] and u w ~ 0.05 [31]. At late time 
(compared to the nucleation time), the bubble profile can be approximated as planar and 
physical quantities depend only on the coordinate z = (r — v w t), the distance to the wall. 
Motivated by realistic bubble wall profiles pj, we take 

. . Vo ( , , 2z \ « / , , 2z \ 

f (z) = — I 1 — tanh — I , a[z) = — I 1 — tanh — I . (9) 
2 \ L w J 2 \ L w J 

The z < (z > 0) region corresponds to the (un)broken phase. Additionally, for simplicity 
we take constant diagonal elements m 2 L R . 



III. QUANTUM BOLTZMANN EQUATIONS 

A. Kinetic Theory Derivation 

Kinetic theory, described by Boltzmann equations, is a useful tool to describe the dy- 
namics of an ensemble of quantum states [28Tl3~T] . For a single species, characterized by the 
distribution function /(k, x) of states with momentum = (wk, k) and spacetime coordinate 
x M = (t,x), the usual Boltzmann equation is (in flat spacetime) 

(d t + v • V x + F • V k ) f(k, x) = tf(k, x) . (10) 

Here, v = k/cut is the velocity and F(k, x) is the force associated with the variation of an 
external potential over length scale L ext . The collision term ^(k, x) characterizes scattering 
interactions, with mean free path L m f p . The Boltzmann picture is valid only in the semi- 
classical limit, where L m f p , L ext ^> L int , where L int = |k| _1 is the "intrinsic" de Broglie 
wavelength 1 . 

In the two-flavor case, a new effect can arise: flavor oscillations due to quantum coherence 
between different mass eigenstates. The relevance of flavor oscillations depends on the 
oscillation length scale L osc oc l/(wi — 002) ~ |k|/Am 2 , where Am 2 = m\ — m\. Clearly, if 



we want to include quantum coherence in our dynamics, Eq. ( 10 ) must be generalized. 



In the context of EWBG, the relevant length scales are given as follows: 



1 This statement follows from the uncertainty principle: Ak Ax > 1. By describing the system in terms of 
a distribution /, it is assumed that states have well-defined momenta, such that Ak <S |k|. At the same 
time, it is assumed that the force term acts locally at x, while the collision term is formulated in terms of 
localized, single scattering interactions. Both assumptions require that states are sufficiently localized in 
position with respect to the relevant scales: Ax <C £ m fp, L cxt - Thus, L m f p , L cxt ^> Ax > 1/Afc ^> Lint- 
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• The typical de Broglie wavelength is L int = |k| 1 ~ T 1 . 

• The external length scale L ex t is set by the wall thickness L w . Previous studies have 
found L w ~ (20 - 30)T" 1 in the MSSM [8] and L w ~ (2 - 40)T~ 1 in extensions of the 
MSSM [331 125] • 

• The oscillation length L osc is determined by the $ mass spectrum. In the thick wall 
regime, CP asymmetries are maximized for L w ~ L osc (see Ref. [16] and the discussion 
in Sec. IV); thus, the L osc 3> L in t case is the most interesting for EWBG. 

• The mean free path satisfies L mfp 3> L int if the couplings yh,R are perturbative. 
Therefore, we assume in our analysis that the following ratios are small parameters: 

Lint _ Li n t Li n t 

^wall = ~~ j ; ^coll = ~p j ^osc = ~p j I / 

l-'w -^mfp J-'osc 

collectively denoted as e. 

In the e<l limit, a Boltzmann-like description of a multi-flavor system is still possible, 
despite the inherently quantum nature of the coherence between states. In this case, one 
must promote /(k, x) to a density matrix: the diagonal elements denote occupation numbers 
of states, while the off-diagonal elements describe coherence between those states. The 
Boltzmann equation becomes a matrix equation: 

(d t + v • V x + F • V k ) f(k, x) = -i [ft k> /(k, x)\ + if [/, /] (k, x) . (12) 



The general structure is nearly identical to Eq. (10), except for two important differences: 

• The free Hamiltonian is now a matrix, f2k = A/|k| 2 + M 2 (x), and gives rise to the new 
commutator term [Q^, /]. 



The collision term /], evaluated explicitly in Appendix|Aj has a non-trivial matrix 
structure involving y and the density matrices for particles (/) and antiparticles (/). 
(Our notation &[f, f] indicates that ^ is a functional of /, /.) 



Although Eq. (12) is covariant under flavor rotations, it is most convenient to work in the 
mass basis, denoted by the subscript m. Rotating to this basis, the density matrix and free 
Hamiltonian transform as 



/(k, x) -> / m (k, x) = u\x)f(k, x)u(x) , n k -> Wk = f ^ k ^ ^ 



(13) 



where = A/|k| 2 + mf(x). The Boltzmann equation becomes 

{dt + v • V x + F • V k ) / m (k, x) = - [ ?Uk + S° + v • S, / m (k, x)] + ^ m [f m , f m ] (k, x) . (14) 

A similar equation governs the evolution of the antiparticle density matrix f m (k,x). The 
quantum Boltzmann equations for f m and f m are derived more rigorously below. The final 
result for our Boltzmann equations is given by Eq. (33). 
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B. Field Theory Derivation 



We derive the Boltzmann equations using non-equilibrium quantum field theory in the 
real time Closed Time Path (CTP) formalism [T7j. The arguments presented here are similar 
to those in our previous work, where we derived the multi-flavor Boltzmann equations for 
purely time-dependent scalar systems |16| . to which we refer the reader for greater detail. 

In the CTP formalism, the basic building blocks are the non-equilibrium Green's func- 
tions, defined here for mass-basis fields 4>i 2 : 



G>(y,y>) = U l (y)ri(y': 



G<(y,y') = U\(y')My) 



G\ J {y,y') = (YUy)<P){y l : 



(15a) 
(15b) 



where T (T ) denotes (anti)time-ordering. These Green's functions obey Schwinger-Dyson 
equations, which along with the free field equations following from the Lagrangian Eq. (J6]) 
imply the equations of motion 



d 2 + m 2 (y) + 2E M (y)fl£ + E M E%) + d^(y)\ G*{y, y' 
GHV, y') \d} + mV) - + Z^(y') - d^(y') 



-i I d 4 z 
-i \d 4 z 



U(y,z)G(z,y') 
G{y,z)n(z,y') 



where the tildes denote matrices in CTP space, 

Qt —Q< 

g = I q> Qt 1 , n 



it -n< 
n> -n* 



(16) 



(17) 



and U(y,y') is the matrix of self-energies that appears in the Schwinger-Dyson equation. 

From these equations of motion, one derives the so-called constraint and kinetic equations 
for the Wigner-transformed Green's functions 



G<(k, x) = / d 4 r e ik - r G} 3 (x + r/2, x - r/2) 



(18) 



where x = (y + y')/2 and r = y — y'. (We will let the arguments distinguish the position 
space Green's function and its Wigner transform.) 
The constraint equation, given by 



2k 2 - ^ ] G*{k,x) 



( { m 2 (x) - 2i k ■ E(x) + S(x) 2 , G*{k, x) } (19) 
+ i{U h {k,x),G^{k,x)} +i{U^{k,x),G h {k,x)} 
+ % - [U>(k,x),G<(k,x)] + % - [G>(k,x),U<(k,x)] \ , 



2 Although our formalism is covariant under flavor rotations, we work in the mass basis since it makes our 
e expansion the most transparent. 
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determines the shell structure of the excitations. The kinetic equation, given by 

2k-d x G*(k,x) = e- i0 ^-2[m 2 (x)-2iA;-£(x) + £(x) 2 , G*(k,x)] (20) 

+ [U h {k,x),G^{k,x)} + [U^(k,x),G h (k,x)] 

+ \ {U>(k,x),G<(k,x)} - \ {U<(k,x),G>(k,x)} S j , 

governs the dynamics of the system 3 . The diamond operator O is defined by 
>v/ r ./, nt,/, ^ 1 { dA dB OA dB \ 

The II functions are now the Wigner-transformed self-energies and G h = (G* — G*)/2. 

Working at a fixed order in a perturbative expansion in couplings Ul,r-, one can express the 
self-energies U(k,x) as functionals of G(k,x). Eqs. (19) and (20) then describe the quantum 
evolution of the Wigner functions G(k,x). However, Eqs. (19) and (20) are formidable to 
solve in practice. Therefore, we simplify them by working to leading non-trivial order in e: 
0(e) in the kinetic equation and O(e ) in the constraint equation. This follows the spirit 
of the effective kinetic theory developed in Refs. [28H3T] . Our power counting in e proceeds 
according to the following rules: 

• Each derivative d x acting on U(x) or m 2 (x) carries one power of e wa n; e.g., £ M is 
0(e wa u). 

• Each factor of the self-energy II carries one power of e co u. This is equivalent to an 
expansion in coupling constants y^R. 

• Each Am 2 carries one power of e osc . In particular, the commutator [m 2 , G^(k,x)} and 
G h are both proportional to Am 2 and are 0(e osc ) 



According to these rules, all terms on the right side of Eq. (20) are at least linear in e. 
Therefore, the kinetic equation implies one more rule: 

• Each derivative d x acting on G^(k,x) carries one power of e. 

The dimensionality of these quantities (e.g., Am 2 , d x U, etc.) is compensated by powers of 
or |k|, taken to be 0(T), to form dimensionless ratios e. Our e expansion, therefore, 
breaks down for infrared modes |k| <T. We neglect this complication since the density of 
states for these modes, k 2 f(k), is suppressed compared to typical thermal modes |k| ~ T. 
Using the above rules, the constraint equation at O(e ) becomes trivial 4 : 

(k 2 -m 2 {x))G^{k,x) = , (22) 



3 Eqs. ^ and ([20| correct typos in Eqs. (25) and (26) in Ref. [IB] . 

4 In general there exists a more complicated shell structure that deserves mention [36] . For free fields 
(setting S,n — >• 0), it is straight-forward to solve the constraint equation to all orders in e osc . In the 
rest frame of the wall, for a given component Gfj{k, x), there exist not two but four shells, which can be 
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where m 2 = (m 2 + m|)/2. Therefore, G^(k,x) must vanish unless 



k° = ±Uk(x) = ±A/|k| 2 + nr(.r) 



(23) 



The two shells correspond to particle (k° > 0) and antiparticle (k° < 0) modes, and the 
two-point functions can be expressed in terms of particle (f m ) and antiparticle (f m ) mass 
basis density matrices as follows: 



(24) 



G > (k,x) = 27i5(k 2 -m 2 ) [8(k°)(I + f m (Kx)) + d(-k )f m (-k,x)} , 
G<(k,x) = 2n5(k 2 -m 2 ) [6(k )f m (Kx) + 9(-k°)(I + f m (-Kx)] , 

where I is the 2x2 identity matrix. 

The Boltzmann equation is obtained from the kinetic equation (20). Working to 0(e) 
we have 



2k-d x G^(k,x) = -[ 



ira (x 



+ 2k- E(x) -U h (k,x), G*{k,x)] 



(25) 



+ i {^m 2 , ^,G<(A;, a ;)} + ^ {n>(fc, x), G<(k, x)} - l - {U<(k,x),G > (k,x)} 



By taking the positive (negative) frequency integrals of G^(k,x), we can project out the 
particle (antiparticle) density matrices: 



/°° dk° - f 

— 2k°G < (k,x), / m (-k,x) = J 





27 



2fc°)G>(fc,x) . 



(26) 



Taking the positive frequency integral of Eq. (25), we arrive at the Boltzmann equations for 
f m . Here, a useful relation is 



2^ 



G < (A;,x) 



°°dk° {2k° 



2n \ 2u\t 



G < (k,x) + 0(e) = t%^- + 0(e) 



(27) 



using Eq. ( |23[ ), according to which the factor (fc°/a>k) is equal to unity (restricted to k° >0), 
modulo 0(e) corrections. These corrections can be neglected since we are working to linear 
order in e and every term in Eq. (25) is already 0(e). 



expressed as (using over-bar for quantities in the wall rest-frame) 

k z = ±(k z i + k Z j)/2 , k z = i(fczj — k Z j)/2 , 



where k zi = \/ k^ — k% — — mf. The k z = ±(k Z i + k Z j)/2 shells describe coherence of states moving 



in the same direction (for i ^ j these are different eigenstates) and reduce to Eq. (23 1 for e osc = 0. The 
other shells k z — ±(k Z i — k Z j)/2, named "non-local coherence shells," correspond to coherence between 
states of opposite momentum (and for i =/= j different mass eigenstates), and are interpreted in terms 
of quantum mechanical reflection |36] . By performing a mode expansion of in terms of free-field 
creation and annihilation operators, it is possible to show that the coherence shells arise from non-zero 
expectation values {aibj} and (at&t); such an effect is also known as zitterbewegung [37] . In the thick wall 
regime (L w >• L- mt ) of EWBG, we expect on physical grounds that occupation numbers associated with 
reflection should be suppressed and we neglect these shells in our analysis. However, their importance in 
thin wall regime (L w < L; nt ) has been emphasized in Refs. [55] , 
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We now evaluate the various terms in the Boltzmann equation. The left side of Eq. (25) 

is 

roc juO 

/ -2fc.a i G<(fc, I ) = (9 t + vV x )/ m (k,x) + 0(e 2 ), (28) 
Jo 2tt 

with velocity v = k/(X>k- The oscillation term is 

roc m,0 

J 27 [ m2(x) ' G< ^ x )] = i^,UKx)} + 0(e 2 ) , (29) 
using the fact that (raf — m 2 ,)/(2Cd k ) = (wn — u; 2 k)- The CP-violating source term is 

roc juO 

/ — [2k-V(x),G<(k,x)] = [E°(x)+vS(x), / m (k,x)]+0(e 2 ). (30) 
Jo 27r 



and the force term is 

rD0 rfA; 1 



/ ^^{^m 2 ,^G<(A; ) x)} = -F-V k / m (k,x) + 0(e 2 ), (31) 
J 2tt 2 



with force F = — V x Wk- (The d^o contribution to Eq. (31) is a total derivative and vanishes 
at the boundaries.) 



The remaining terms in Eq. (25), arising from the self-energy IT, give two important con- 
tributions (see Appendix IaJ. First, the [H h , term yields a medium-dependent, forward- 
scattering correction to the mass matrix. For the interaction given in Eq. (|2]), assum- 
ing the A bosons are in thermal equilibrium, this correction gives the thermal mass shift 
m 2 L R — > m 2 L R + Ul,rT 2 /24. This shift can be incorporated directly into m(x) and S M (x). 
The remaining collision term 

/°°dk 1 
— - ({ll>(k,x),G<(k,x)} - {ll<(k,x),G>(k,x)}) (32) 

corresponds to scattering (<pA -h- <pA) and annihilation (00^ -h- AA) processes in the plasma. 
To summarize, the quantum Boltzmann equations are 5 

(u- d x + F ■ Vk) /m(k, x) = - [iu) k + u ■ S, f m (k,x)] + ^ m [/ m , f m ](k,x) (33a) 
(u-d x + F- V k )/ m (k,aO = +[iw k -«-E, / m (k,a;)] +# m [/ m ,/ m ](k,x) (33b) 

for the (anti)particle density matrix f m (f m ), with = (l,v) and d% = (d t , V x )- 

These equations are identical in structure to the usual single-flavor Boltzmann equations, 
with two additional ingredients. First, the term [iu^, f m ] gives rise to ($£,$#) flavor oscil- 
lations. Second, the [u ■ E, f m ] term is the CP-violating source, due to spacetime-dependent 
mixing. This term is a "source" because it does not vanish when f m , f m are in equilibrium; 
furthermore, it violates C and CP symmetries, as we show below. Lastly, it is straightfor- 



ward to show that Eqs. (33) are consistent with the continuity equation for the total <3>£ + $^ 



charge, Tr / d 3 k/(2ii) 3 (u ■ d x )(f m - f m ) = 



5 We follow a convention where the antiparticle density matrix obeys the same flavor transformation rule 



/ — > U> f m ll as / in Eq. ( 13 ). If we evaluated the fields </>, $ in terms of creation/annihilation operators, 
the density matrices would be (f m )ij ~ (a] a i) and (fm)ij ~ (tybj)- The swapping of i,j between f m 



and f m is the reason for the sign flip between the terms [wk, /m] and [uJk, f m ] in Eq. (33 1. If we chose 
the alternate transformation convention / — > U T f m U* , we would have (f m )ij ~ (bjbi), no sign-flipped 
[uk,/ m ] term, and S p replaced by £ M *. 
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C. C and CP violation 

It is insightful to consider how C and CP violation are manifested in the Boltzmann 



equations (33). (See Ref. [IB] for a complementary discussion at the Lagrangian level.) 
Under C, the density matrices transform as 6 

f m (k,x) e i y^(k,x)e- i \ (34) 

where T] = diag(r/i, 772) are arbitrary phases. The Boltzmann equations are C-symmetric if 

e i„ s t e -i V = _ s ^ ( 35 ) 

for some choice of 77. Therefore, from Eq. ([7]), C violation requires sm9 7^ and d^a 7^ 0. 
Under P, the density matrices transform as 

/ m (k,x) e^7m(-k,x')e-^, (36) 

where fj = diag(f/!, f/ 2 ) are again arbitrary phases, and x ,fl = (t, — x) is the P-inverted 
coordinate. In general, inhomogeneous background fields break spatial symmetries, such as 
P. However, a spherical bubble centered at r = is invariant under P: therefore, we have 

m 2 ( x ) = m 2 ( x ') ; U(x) = U(x') , F(x) = -F(x') , £"(x) = (S°(x'), -E(a/)) , (37) 

using the fact that V x = — V x '- The Boltzmann equation for f m transforms under P into 

(0 t -vV x ,-F- V k )/ m (-k,x') 

= - [^(xO + EVO-v-S^'), / m (-k,x')] +^ m [/ m ,/ m ](-k,x') , (38) 



taking 771,2 = 0. Therefore, setting k — > — k and relabeling x' — >■ x, we find that Eqs. (33) are 
invariant under P. 

In summary, the Sakharov conditions of C and CP violation are realized if the bubble 
wall induces flavor mixing (sin^ 7^ 0) and a spacetime-dependent phase (d^a 7^ 0) in the 
two-scalar system. For a spherical bubble, C and CP violation are equivalent, since P is 
conserved. An aspherical bubble will in general violate P, and therefore CP, but clearly this 
is insufficient by itself for EWBG if C is conserved. 7 



IV. SOLUTION TO THE FLAVORED BOLTZMANN EQUATIONS 



In this section, we solve the flavored Boltzmann equations numerically, organizing our 
discussion as follows. In Section IV A we apply the formalism of Sec. Ill to our EWBG 
toy model of mixing scalars $ in a bubble wall geometry, and we describe what are the 



quantities of interest for EWBG. In Section IV B we show that EWBG shares a physical 



analogy with spin precession in a varying magnetic field. In Section |IV C| we describe our 
numerical methods for solving the Boltzmann equations. We present our numerical results 



and discuss the implications for EWBG in Section IV D Our main conclusions are: 



The reason for the transpose is our convention for f m given in footnote [Ej In the alternate convention, 



one would have /„ 



7 Later we will take the planar limit of the bubble wall as given by Eq. which apparently violates P. 
However, Eq. (J9J) really only describes one "side" of the bubble in the planar limit, with the other side 
infinitely far away, hiding its true parity invariance. Under P, Eq. ([§ remains unchanged. 
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• We demonstrate what are the key charge transport dynamics of EWBG: CP- violating 
flavor mixing and coherent oscillations in a spacetime-dependent background, colli- 
sional damping that destroys this coherence, and diffusion of charge into the unbroken 
phase. Our solutions are exact, without any ansatz for the functional form of the $ 
density matrices /, /. 

• We find a resonant enhancement of CP- violating charge generation for mi ~ tur. 
Similar resonances were discussed previously in Refs. p2H2l], but our present work 
establishes its true origin on more theoretically sound footing. The width and height 
of this resonance are controlled by the ratio L osc / L w for typical thermal modes. 

• CP-violating charge does diffuse into the unbroken phase. This can provide a po- 
tentially large enhancement of charge production compared to previous treatments 
of EWBG flavor oscillations in Ref. [Sj, which found that diffusion of the oscillating 
species was quenched. We identify the reason for this discrepancy in Sec. [V] finding 
it is due to the breakdown of the power counting expansion followed in Ref. [H] . 

To clarify our strategy for power counting solutions to the Boltzmann equations, we dis- 
tinguish between power counting terms in the equations themselves, and explicitly expanding 
the solutions for G^(k, x) perturbatively in e. We perform the former but not the latter. For 
the equations of motion at a given order in e, we will solve for G(k, x) exactly as a function 
of e. We define the 0(e n ) solution for G(k,x) as the exact solution of the C(e n ) equation of 
motion. To be precise, we will obtain the exact solutions of the constraint equation at O(e ) 
and the kinetic equation to 0(e), an appropriate strategy to solve for the leading nontrivial 
deviations of the distribution functions away from equilibrium. This method avoids making 
any a priori ansatz about the functional form or power counting of G(k,x). 8 



A. EWBG setup 

For EWBG in a late time regime (compared to the bubble nucleation time), significant 



simplifications arise in solving the quantum Boltzmann equations (33). First, we neglect 
the wall curvature, treating the spherical bubble as a planar wall, where z < corresponds 
to the broken phase inside the bubble, as given by Eq. Second, we look for steady 
state solutions in the rest frame of the moving wall (with v w <C 1) for the $ particle and 
antiparticle density matrices /(k, z) and / (k, z), assumed to be a function only of z = r—v w t, 
the coordinate normal to the wall. Additionally, the force term F vanishes in our model 



since ml + m| is constant. Thus, Eq. (33) becomes 9 

v rcl d z f(k,z) = -[iu k + v iel E, f(k,z)] +tf[f,f](k,z) (39a) 
v iel dj(k, z) = [i u k - v rcl S, /(k, z)] + if [/, /] (k, z) , (39b) 

with 

tfea(k) = ^-v w , u; k (z) = ( Wlk Q (z) ^ ) , E(z) = U\z)d z U{z) , (40) 

8 Furthermore, it is only by this method that one obtains solutions that equilibrate properly at late time 
(or far from the bubble wall). Otherwise we run into the problem of "secular terms" that grow large and 
spoil equilibration at late time [52]. Our procedure will not encounter any such problematic terms. 

9 Henceforth, we work in the mass basis and drop the m subscripts whenever it does not lead to ambiguitues. 
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where v re \ is the velocity with respect to the wall and n is the unit vector normal to the 
wall. From now on, the collision term 'rf includes a factor of g* to model the true number 



of degrees of freedom in the electroweak plasma and is given by Eq. (A4). Given a set of 



input model parameters, we want to solve Eq. (39) subject to the boundary condition that 
far from the wall the solutions reach equilibrium: 

lim /(k,z),/(k,z) = lim r(Kz) } n(k,z)=( nB{uJ f z)) ° ) . (41) 



Since Eq. (39) has azimuthal symmetry with respect to n, the density matrices /, / depend 
only on the momentum variables k = |k| and cos$k = k • n, and are independent of the 
azimuthal k angle. 



After solving Eq. (39), we compute total charge asymmetries that are directly relevant 
for EWBG. The charge current matrix (in the mass basis) can be defined equivalently in 
terms of normal-ordered fields, Green's functions, or density matrices: 

Z(-r)u = / -^k* (Gg(*. x) + G>(k, x)) 

d 3 k k» 



(27TJ 

SwFlL {UKx)-UKx)) io . (42) 



The analogous flavor-basis charge current is j£ = U U\ and this coincides with in the 
unbroken phase (z>0). The quantities 

ni(z) = {f m )n , n 2 (z) = {j° m ) 2 2 , n L (z) = (j$) LL , n R (z) = (j^rr (43) 

are the spacetime-dependent total charge densities of <pi, 4>2, and states, respectively. 
In the unbroken phase, = til,r- In addition, we define the integrals 

POO 

I C l P r= / dzn L>R (z) (44) 

as the total charge in the unbroken phase. This is a useful global measure of CP violation 
and a suitable proxy for the baryon asymmetry itself. For example, in squark-driven EWBG 
scenarios [53], where $ = (t L ,tji), and will be converted into left-handed quark 
charge through fast gaugino- and Higgsino-mediated processes, respectively, thereby sourcing 
baryon generation through electroweak sphalerons (which are only active in the unbroken 
phase). Therefore, baryon number will be directly proportional to Ilr- (A more precise 
statement requires a "realistic" model such as [33], beyond the scope of this work.) 



B. Magnetic analogy 

At this point, it insightful to introduce an analogy with spin precession in a varying 
magnetic field [38]. The 2x2 Hermitian density matrices /, / can be expressed by Bloch 
decomposition as four- vectors p = (po, p) and p = (po, p): 



f(k,z) = I Po (k, z) + a ■ p(k, z) , f(k,z) = I Po (k,z) + <t -p(k,z) , (45) 
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where er = (a±, 02, 03) are Pauli matrices and I is the identity. po(k, z) represents the total 
occupation number of all 0i + 02 particles, while the "polarization" vector p(k, z) describes 
the density matrix for the internal flavor degrees of freedom, for a given momentum k (and 
similarly with p(k, z) for anti-particles) 10 . One can also decompose the collision term as 

— <€ [f, f] (k, z) = - 1 D [p, p] (k, z) - a- ■ D [p, p] (k, z) (46) 

^rel 

to define the damping vector D = (_D ,D). 

In the polarization vector language, Eq. (J39j) becomes 

d z p(k,z) = (B + B s )xp(k,2)-D[p,i)](k, Z ) (47a) 
d z p(k,z) = -(B -B E ) xp(k,z)-D[p,p](k,z) (47b) 

with effective magnetic field given by 

(48a) 

2 sin 2 9 a') ,(48b) 

where the primes ' denote derivatives d/dz. The equations of motion are strongly suggestive 
of spin precession in a magnetic field. Here, flavor polarizations p, p play the role of spin 
and precess around an effective magnetic field (Bo ± Be) in flavor space. This describes 
coherent flavor oscillations. The collision term D destroys coherent oscillations by damping 
the precession. The total particle/antiparticle occupation numbers obey Boltzmann-type 
equations with collisions 

d z p Q (k,z) = -D Q \p,p](k,z) , d z p (k,z) = -D [p,p](k, z) . (49) 

This analogy provides an intuitive framework to understand the behavior of our numerical 
results. We illustrate this picture in Fig. [TJ The qualitative features of the solutions are 
controlled by two ratios of scales: 

• Oscillation vs. wall length: the ratio L osc /L w , with oscillation length L osc = 
2tt v ie i/(u)xk — ^2k) and wall thickness L w (which determines the z-dependence of the 
effective magnetic field), controls how quickly the magnetic field varies on the intrinsic 
time scale of the system. In the adiabatic regime (L osc /L w 1), the polarization 
vector tracks the magnetic field with small amplitude precession. In the non-adiabatic 
regime (L osc /L w > 1), the polarization vector precesses with large amplitude. 

• Collisional mean free path vs. wall length: the ratio L co \\/L w , with typical 
collision mean free path L co \\, controls how fast the CP asymmetry is damped away and 
equilibrium is restored away from the wall. In the overdamped regime, (L co \\/L w 1) 
precession is efficiently damped away, allowing no CP asymmetry to develop. In our 
analysis, we consider L co \\/L w > 1; interactions do not affect sizably the sourcing of 
CP asymmetries, but re-establish equilibrium away from the wall (\z\ L w ). 



We denote unit vectors in the internal flavor space by (x, y, z), not to be confused the spatial coordinate 
z with respect to the wall. 



B (k,z) 
BeW 



0,0 



UJvAz) 



^rcl(k) 



(2 sin a 9' + sin 29 cos a a', —2 cos a 9' + sin 29 sin a a' , 
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FIG. 1: Precession of flavor polarization vectors about effective magnetic fields in the mass basis. In 
equilibrium, p = p = (0, 0, rig^i)— ng(w2))/2 and the magnetic fieWBo = (0,0, (u)%— 0J2)/v r ei) both 
point along the z direction, and there is no precession. When the bubble wall turns on, it induces an 
additional magnetic field B^, causing p, p to precess around Bo ± B^, which corresponds to flavor 
oscillations. If a' 7^ so that there is CP-violation, p, p develop different z components, generating 
CP-asymmetric diagonal densities. In the adiabatic regime, Bx; <C Bo, the angle of precession is 
small, the system remains near equilibrium, and large CP asymmetries are not generated. In the 
non-adiabatic regime, B^ > Bo, the precession angle is large, and large deviations from equilibrium 
and large CP asymmetries can arise. The precession is damped by collisions with A bosons in the 
thermal bath, leading the system back to equilibrium at late time. 



These above concepts will serve as organizing principles for the discussion in the following 
sections. 

Lastly, we note 

1 f d 3 k 

nij(z) = - j -j—^ (po(k,z) ±p z (k,z) -p (k,z) =|=p z (k,z)) . (50) 

That is, the charge densities are determined by the differences p z —p z and (to a lesser extent, 
as it turns out) p — p . In equilibrium, p and p are aligned with B , and p Q = p = 
(n B (uJik) + n B (u 2 ] i ))/2 and p z = p z = (n B (u; lk ) - n B (w 2k ))/2. 
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L w = § v w = ± v = T' a = f m L = 2.2T m R = 2T y L = 1 y R ■ = 0.75 5 * = 200 



TABLE I: Parameters that define the "baseline model" used to illustrate the main features of the 
numerical solution. All dimensionful parameters are expressed in units of the temperature T or its 
inverse. In addition, we assume thermal masses (itia/T) 2 = (y^ + y R )/12 for the scalar field A. 

C. Numerical approach 



The Boltzmann equations, equivalently described by Eq. (39) and Eq. (47), are a system 
of 8 coupled integro-differential equations, due to the collision term coupling together modes 
of different momenta. In order to make this problem tractable, we discretize k = |k| and 
cos^fc into Nk and N& bins, within the ranges 

< k < k max , -l<cos0 fc <l, (51) 

evaluating the discretized (k, cos$k) at the cental value of each bin. After binning, we have 
a system of 8 x N k x N# coupled first order ordinary differential equations with boundary 
conditions. We solve this system of equations using the "relaxation method" [39] . 



We impose the boundary conditions (41) as follows: for right-moving modes (v Te \ > 0), 
we set /(k, z_) = / cq (k, z_) with z_ < far in the broken phase, and for left-moving modes 
(f re i < 0), we set /(k, z + ) = / cq (k, z + ) with z + > far in the unbroken phase. These split 
boundary conditions are required on physical grounds: the collision term equilibrates the 
density matrix in the positive time direction, which for right (left)-going modes is the positive 
(negative) z direction. So we have only to impose equilibration as a boundary condition in 
the negative time direction, i.e. negative (positive) z for right-(left-)moving modes. The fact 
that right (left)-going modes equilibrate again at late time, for z — > z + (z — > zJ), provides a 
non-trivial check on our numerics. 

Within the baseline model parameters in Table [TJ we have performed a number of stability 
checks against different choices of /c max , Nk, N#, z±. We required that quantities of physical 
relevance, e.g. n LjR (z), remain stable at the percent level. We find that \z±\ = 5000/T and 
fc max = 8T are acceptable values. 11 Moreover, binning as coarse as N k = 4, N$ = 6 produces 
stable charge density profiles. Our results below have N^ = N$ = 8. 

Given our discretized solutions for f(k,z) and f(k,z), we compute the charge currents 



and densities in Eq. (42) by converting the continuous integrals into discretized sums in the 
usual way: / d 3 k/(2nf-^ k ma J (2n 2 N k N$) E fe , C os^ fc2 - 

D. Numerical results 

1. Distribution functions 

Numerical results within the baseline model defined in Table |T] are presented in Fig. [2j 
The figures show the z dependence of Po,x,y,z (left panels) and and po, x , y ,z (right panels). The 



ii 



Decreasing yL,R and/or increases i C oii, thus requiring larger values of \z±\. 
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FIG. 2: Numerical results for particle and anti-particle density matrix for [top] a typical adiabatic 
bin: k/T = 0.5,cos??£; = 0.875, and [bottom] a typical non-adiabatic bin, k/T = 3.5, cos = 
0.875. Left panels display particle density matrix in Bloch decomposition: po(k, cos z) 
and p X} y^ z (k, cosi?^, z). Right panels display anti-particle density matrix: po(k, cosi?/-, z) and 
Px,y,z(k, cos i?^, z). Solid lines represent full numerical solutions, while dotted lines represent lo- 
cal thermal equilibrium results for diagonal components. See text for additional details. 
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top four panels correspond to a typical adiabatic bin with k/T = 0.5, cos$& = 0.875 and 
L OS c( z — 0)/L w ~ 0.16, while the bottom four panels correspond to a typical non-adiabatic 
bin with k/T = 3.5,cos-# fc = 0.875 and L osc (z = 0)/L w ~ 1.37. 12 In all plots the solid 
lines represent the full numerical solutions, while the dotted lines represent local thermal 
equilibrium results. 

As anticipated, the qualitative behavior of Fig. [2] can be readily understood through the 
magnetic analogy. Let us neglect for a moment the effect of interactions with the thermal 
bath. In the collisionless limit different momentum bins are decoupled: po and po do not 



evolve, while p and p precess about effective magnetic fields as per Eqs. (47). Equilibrium 
boundary conditions for left- and right-moving modes imply that for \z\ ^> L w both p and 
p point along the z axis in flavor space. This is a stable configuration as long as Be = 0. In 
proximity of the phase boundary, the non-vanishing Be tends to push p and p out of their 
stationary state, triggering the precession around the ^-dependent fields B ± B s . 

In the adiabatic regime (L w ^> L osc ), the polarization vectors p and p effectively track the 
magnetic fields Bo±Be (with a small precession amplitude that vanishes in the L osc /L w — > 
limit). As a consequence the solution tracks very closely the local thermal equilibrium. On 
the other hand, in the non-adiabatic regime (L w < L osc ), when the magnetic field changes 
on length scales comparable to or smaller than the oscillation scale, the polarization vector 
lags behind the magnetic field and begins precessing with a large amplitude. (In absence of 
collisions the precession persists away from the phase boundary, \z\ ^> L w .) The amplitude 
of oscillations increases with L osc /L w , as is evident from Fig. |2j in the non-adiabatic regime 
the system is pushed out of equilibrium more efficiently by the passage of the bubble wall. 
Collisions and pair processes play an essential role in relaxing the density matrices back to 
equilibrium away from the phase boundary, as evident from the plots in Fig. [2j 

In the non-adiabatic regime, CP-violating effects show up more prominently in the evo- 
lution of the density matrices. In the CP-conserving limit (a'(z) = 0), the effective mag- 
netic fields Bo ± Be are confined to a plane defined by z and (sin ax — cos ay). The 
evolution obeys the CP invariance condition f(k, cos d^, z) = e tr, f T (k, cosfik, z)e~' tv , where 
i] = diag(a, — a)/2, or 

Px - ip y = e~ ia (p x + ipy) , pz=Pz, Po = Po ■ (52) 

In presence of CP violation (a'(z) ^ 0), Bo ± Be are not confined to this plane, so that the 
dynamical evolution leads to an angle between p and Bq + Be different from that between 



p and Bo — Be- This leads to a violation of the conditions (52) and generation of flavor- 
diagonal CP asymmetries p z — p z ^ (and eventually, through collisions, p — p ^ 0). 
The CP asymmetries vanish in two limits: (i) L osc /L w 1, because during the resulting 
adiabatic evolution the polarizations track closely the magnetic fields and so end up in CP- 
symmetric thermal equilibrium, and (ii) L osc / L w 3> 1, because, as discussed in Ref-PS], then 
the magnetic field varies so fast that precession becomes sensitive only to the initial and final 
values of Bo ± Be, which define a plane. CP asymmetries are maximal for L osc /L w ~ 0(1). 

In Fig. [3j we illustrate the "anatomy" of how CP asymmetries are sourced for different 
momentum bins. Since L osc ~ 2iik / \m\—m\\, we expect greater asymmetries to be generated 



12 In general the adiabaticity L osc /L w is controlled by fc,cos^fe, mi 2 , and L w . Larger values of k and 
| cos$fc| and smaller mass splittings increase L OSCl thus leading to increasingly non-adiabatic evolution for 
fixed L w . Once the underlying model parameters have been fixed, the adiabaticity is controlled only by 
k and cos . 



19 




0. 



010 



0. 



012 



-100 







100 



200 



300 




FIG. 3: Diagonal CP asymmetries k 2 (f\\(z) — fu(z)) (in units of T 2 ), for three representative 
bins. The long-dashed line refers to a typical adiabatic bin (k/T = 0.5, cos = 0.875), the short- 
dashed line to a mildly non-adiabatic bin (k/T = 1.5, cos = 0.875), and the solid line to a typical 
non-adiabatic bin (k/T = 3.5, cos = 0.875). 

for non-adiabatic bins corrsponding to larger values of k. We plot the diagonal <f>\ CP 
asymmetry k 2 (fu — /n), for three representative (k/T, cos$fc) bins: the long-dashed line 
refers to a typical adiabatic bin (0.5, 0.875), the short-dashed line to a mildly non-adiabatic 
bin (1.5, 0.875), and the solid line to a typical non-adiabatic bin (3.5, 0.875). We weight each 
asymmetry by the phase space factor k 2 , so that Fig. [3] represents the contributions of each 
bin to the total charge density ni = J d 3 k(f n — /n)/(27r) 3 . The plots clearly illustrate: 

(i) In the vicinity of the wall (z < L w = 20/T), the largest asymmetry is generated for 
non-adiabatic momentum bins. 

(ii) Collisions establish kinetic equilibrium away from the wall (L w < z < 100/T), by 
redistributing charge among bins. Far from the wall (z > 100/T), the density matrices 
are well-described by equilibrium distribution functions with a non-zero, spacetime- 
dependent chemical potential. 

This picture holds as long as L co u > L w , which is verified in the baseline model. A qualitative 
difference would arise in the case in which L co \\ < L w . In that case we expect a suppression 
of the CP asymmetries [16], because collisions are so frequent that they break the coherent 
evolution needed for a manifestation of CP-violating effects: then flavor oscillations cannot 
play a significant role in generating a CP asymmetry. 
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FIG. 4: Numerical results for the net charge densities n\{z) (solid line) and ri2{z) (dashed line), in 
units of T 3 , within the baseline model specified in Table^ In the unbroken phase (z > L w ) n\ = til 
and 77,2 = n R (mass and flavor basis coincide). This plot illustrates the existence of diffusion tails 
in the unbroken phase. The relative size of the diffusion tails for n\ and n 2 is as expected, given 
that yi > Ur in the baseline model. 



2. Charge densities and currents 

So far we have presented results for the mass-basis density matrices /, / (in the Bloch rep- 
resentation). In applications to EWBG, one is interested in the behavior of flavor-diagonal 
CP-violating charge densities Ul^z) in the unbroken phase. In Fig. [4]we present numerical 
results for the densities n\(z) and n 2 (z). (For z^> L w one has n L = ri\ and n R = n 2 .) The 
plot in Fig. [4] clearly illustrates the existence of diffusion into the unbroken phase. Once 
generated by CP- violating oscillations within the bubble wall, <&l,r charge diffuses into the 
unbroken phase, where mass and flavor eigenstates coincide and flavor oscillations no longer 
occur. The smaller diffusion tail in n\ — Ul compared to n 2 = tir is due to the fact that 
Vl > Ur in the baseline model, and so the mean free path for is shorter. On the other 
hand, in the broken phase where flavor and mass eigenstates do not coincide, flavor-sensitive 
collisions (yi ^ yn) lead to fast flavor equilibration (til — tir — » 0). Since total charge con- 
servation and causality imply + tir — > far from the wall, flavor equilibration has the 
effect of driving both diagonal densities to zero 13 . 

Another interesting dynamical question involves the onset of the diffusion regime. Within 
our model we can evaluate the current densities and check whether there is a regime in which 
they satisfy the diffusion ansatz ju(z) = —Di'Vrii(z), with diffusion constants Z)j. We have 
found that this ansatz is fairly well satisfied for large z far from the wall (z > 300/T) for 
some numerically fitted constant value for Di. Of course, in principle one should calculate 



13 This has essentially the same effect as the T m rates introduced in the diffusion equation treatment of this 
problem [25H27] . although the physical mechanisms are not identical. 
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FIG. 5: Resonant enhancement of , the total integrated $t charge asymmetry in the unbroken 
phase [see Eq. (44)/, as a function of the mass parameter ttil (with fixed m^/T = 2) for different 
values of the wall thickness L w : L w = 10/T (long-dashed curve), L w = 20/T (solid curve), L w = 
40/T (short- dashed curve). All other input parameters are as in Table^ 



the diffusion coefficients from the collision terms. Nevertheless, this observation implies 
that a simplified treatment in terms of diffusion equations for flavor-diagonal densities with 
appropriate "oscillation-induced" sources might lead to satisfactory results. We leave a more 
detailed investigation of this issue, including calculation of the diffusion constants in this 
model, to a forthcoming paper. 



3. Resonant enhancement of CP asymmetry 

Up to this point we have presented results for one particular point in parameter space, 
defined by the baseline model (Table [T]). In phenomenological applications to baryogenesis, 
one would like to identify those regions of parameter space in which the CP asymmetries (and 
eventually the baryon asymmetry) are maximized. To this end, a useful global measure of 
CP violation and a proxy for the baryon asymmetry itself within the toy model is provided 



by the integral J^ p [see Eq. (44)], which represents the total $l charge in the unbroken 
phase. 

From the discussion in the preceding section we have learned that the size of the CP 
asymmetry is controlled by the ratio L osc /L w of oscillation length to wall thickness. Here one 
should really think about L osc as a thermally averaged oscillation length, or the oscillation 
length evaluated at a typical thermal momentum k ~ 3T. For fixed Vq, L osc is controlled by 
the mass splitting — rriR. In Fig. [5j we plot versus m^/T, for fixed mji/T = 2 and 
all other parameters as in Table [Tj The dramatic resonant feature at mi = itir is interpreted 
in terms of non-adiabatic dynamics discussed in the previous section: for mi ~ ttir the 
average oscillation length L osc is maximized, implying that more momentum modes evolve 
non-adiabatically and therefore develop larger CP asymmetries. Fig. [5] demonstrates in a 
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consistent framework of flavor mixing the resonant baryogenesis regime previously discussed 
in the context of perturbative mass-insertions [T8l IT9| 1231 [23] or perturbative insertions of 
mass gradients [20H22], and places the origin of this resonance on a firmer theoretical footing. 

Finally, one can also study the dependence of the resonant enhancement of on other 
model parameters, such as the wall velocity v w , the coupling constant jjl, and the wall 
thickness L w . Decreasing (increasing) v w and yi increases (decreases) the size of the diffusion 
tail in ul(z), and hence leads to a larger (smaller) 1^ '. The dependence on L w is more subtle 
than a simple overall scaling, because it affects both the peak and width of the resonance, as 
illustrated in Fig. [5j The resonance width is determined by the condition L osc /L w > 0(1). 
Numerically we find considerable resonant enhancement for — < 10/ L w (see Fig. [5J. 
Moreover, as L w changes, the number of momentum modes that evolve non-adiabatically 
also changes, thus changing the overall peak of the resonance (decreasing L w leads to larger 
peak value for /£ p ). 

V. COMPARISON WITH PREVIOUS APPROACHES 

In attempting to obtain a more tractable, analytic solution, all previous treatments have 
employed certain approximations to decouple the diagonal and off-diagonal components of 
the kinetic equations. In this section, we provide a detailed comparison of our results with 
the work of Ref . [9] . That work provided the first derivation of the coupled two-flavor kinetic 
equations using the gradient expansion and treatment of flavor oscillations in EWBG. Earlier 
works neglected quantum coherence implicitly by projecting onto diagonal densities within 
a diffusion- type ansatz. 

Our treatment and that of Ref. [9] differ at the stage of power counting and solving these 
equations. The two primary differences are: 

• Diagonal densities: The power counting of Ref. [9] leads one to neglect the diagonal 
components of the source [S, /] so that /n,22 do not depart from equilibrium at first 
order in 1/(L W T). Effectively, this prevents any CP-asymmetry in the flavor-diagonal 
densities ti^r (generated in the bubble wall) from diffusing into the unbroken phase. 
In contrast, in our treatment we find that in a consistent power counting scheme 
deviations of 711,22 from equilibrium are sourced by nonzero /i 2 , and then diffuse deep 
into the unbroken phase, where the mass and flavor bases coincide. 

• Off-diagonal densities: The off-diagonal density /12 in Ref. [9] is sourced only by 
equilibrium diagonal densities, and its approach to equilibrium is described with a 
phenomenological ansatz for the collision term. While in some regimes of parameter 
space these simplifications capture the qualitative behavior of /12 fairly well, quanti- 
tatively they lead to 0(1) deviations from the exact /12. In our treatment, we account 
for all contributions to /12 from the source and collision terms at leading nontrivial 
order in our power counting. 

Ultimately, the approximations of Ref. [9], applied to our toy model, lead one to neglect 
diffusion and result in a substantial underestimation of charge in the unbroken phase, com- 
pared to our exact numerical treatment, as illustrated dramatically below in Fig. [7j The 
reasons for this discrepancy are explained in detail below. 
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A. Source and Collision Terms 



The coupling of the different components of the 2x2 distribution functions /(k, z) arises 
from two sets of terms. In the mass basis, one is the CP-violating source proportional 
to The other is the collision term c €. Our treatment accounts for the full coupled 

structure of both terms, while Ref. [9] argued that the evolution of diagonal and off-diagonal 
densities could be decoupled in the source and collision terms. 

Beginning with the source term, consider the components of its matrix structure: 

(S12/2I — 5] 2 l/l2 —^12 (fll — /22J + (^11 — S22) fl2 \ 

• (53) 
S21 (/n — ^22) — (Sxi — S22) /21 S21/12 — £12/21 / 

In the diagonal entries, we find that the off-diagonal distribution functions /i2,2i source the 
diagonal distributions /n,22 through £12,21- In the off-diagonal entries, /n,22 feed back to 
act as sources for J \ 2 ,2i- In our work, we have not made any a priori assumptions about 
the scalings of the /'s with gradients of the external field (e wa u) and thus accounted for the 
full coupled evolution. In Ref. [9], however, all deviations of /y away from their equilibrium 
values were power counted as 0(e wa n). That is, 

hbM \ . f 5 fu /12 



V n B (u 2 )J V/21 0/22/ 
where 5/11,22 and /i2,2i ~ C(e W aii)- According to this counting, the source term is 

-E12 {n B {ui) - n B (u 2 )) 

£21 ( n BM ~ n B M) 



Ref. [9] 



+ 0(eU) . (55) 



Ref. [D] argued that, working at C?(e wa u), the C(e^ a n) terms could be neglected. 

Here is the crucial point: in the power counting of Ref. [9], there is no CP-violating source 
for the diagonal densities. Furthermore, although there is a source for the off-diagonal 
densities, the different components of /y are decoupled; CP violation in ^12,21 does not 
feed into /n,22- By rotating to the flavor basis, one has ul,r = =Fsin2# / d 3 k Re[(/i2 — 
,/i2)e"*]/(27r) 3 . Charges til,r vanish in the unbroken phase, since 8 = 0. No diffusion exists. 

Now we consider the collision term for /i 2 . The coupled structure of ^ was simplified in 
Ref. [9] by making a simple phenomenological ansatz 

Ref. 0: Ci 2 (k,z) = -T 12 fi 2 (k,z) (56) 

for some constant relaxation rate Ti2, which was estimated to be ~ aT where a is the cou- 
pling strength of the dominant interaction of the species. As we discuss in Appendix |Bj there 
is indeed a part of the full collision term that takes the form of this ansatz, although with 
a k-dependent rate Ti2(k, z). We will consider in the Appendix the full set of contributions 



to the collision term, and the conditions under which the ansatz Eq. (56) may be justified. 
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B. Power Counting of Off-Diagonal Solution 



We now consider solutions for /, / in our toy model, following the procedure of Ref . [H] , 
by making the assumptions described above for the source and collision terms. One obtains 
the decoupled equations for the off-diagonal densities, 



Vrc\d z + i(ui - u 2 ) + f re i(Sii - S 22 ) + T 12 /i 2 (k, z) = v re iS 12 (57a) 
Vreidz - i{oJ\ - u 2 ) + u re i(£n - £ 22 ) + ^12 /i2(k, z) = v re \Si2 ■ (57b) 



where 



S 12 (k,z) = T, 12 (z) (n B (u lk ) - n B (u2k)) 



(58) 



The effective sources for particles and antiparticles coincide, as they only involve the diagonal 
equilibrium densities = 
will adopt a constant value Ti 2 
and cos^i, = l. 14 



For the relaxation rates Ti 2 ,Ti 2 in Eq. ( [57] ), we 

3T 



r 12 = 0.024 T, obtained by evaluating Eq. M) at k 



The simplified Eq. (57) can be solved explicitly, giving [9] 



fi 2 (k;z) = / dz S 12 (k,z) exp < - I dz" 

J —oo L J z' 



+ (En - S 22 ) + ^ 

frel ^rel 



• (59) 



One can notice from this solution, however, that power counting / 12 as C(e wa u) breaks down 
in the non-adiabatic regime. Now, it is true that E 12 which explicitly appears in the integrand 



of Eq. (59) is C(e wa u). However, the integration measure also must be power counted. It 
is determined by the shortest length scale among the wall length L w , the oscillation length 
L osc or the collision length L co u ~ I/T12, since that will determine the effective range of 
integration. Supposing the latter to be very long, we can just compare L w and L osc . In the 
non-adiabatic regime (L w < L osc ) L w is the shortest scale; the effective range of integration 
in Eq. (59) is L w since £ 12 (z') is nonzero only in the bubble wall region. However, we note 
that L w ~ l/e wa n and f\ 2 in Eq. (59) should be power counted as 



/12 ~ (range of integration) x S 



1 



12 



^wall 



x Cwall 



c 

"wall ' 



(60) 
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invalidating the estimate in Ref. [9] that f\ 2 ~ e wa n. J 

In the adiabatic regime, the power counting of Ref. (9] is formally consistent, as explained 
in Appendix [C] As a result, Eq. (59), sourced only by the equilibrium diagonal distributions 



in Eq. (58 ), is actually a fairly good approximation for the true solution for /12, as illustrated 



below in the top two panels of Fig. m 



14 



In general T^k, z) in Eq. (B4| varies considerably with k (it decreases by a factor of four over the range 



< k/T < 9) but has a weak dependence on and z. 
15 In our power counting scheme, S12 in Eq. (58 1 is also order e OS c, so one may ask if the terms with /i2,2i 
in the source Eq. ( 53 1 may be counted as order e 2 and be dropped, mimicking Ref. [9] . However, as we 
show in Appendix O if one power counts this way then every term in the kinetic equation is order e 2 , 



and every term in Eq. (53) must be kept, thus keeping the equations coupled. 
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The above observations can be easily understood from Fig. [TJ which illustrated the mag- 
netic analogy for precession of flavor polarization vectors. In equilibrium, |p| starts out 
proportional to n B {ui) — n B (u 2 ) ~ e osc . Also the initial magnetic field B ~ oj\ — u 2 ~ e osc . 
When the wall B s ~ e wa n turns on, the precession angle is thus of order 6 B ~ e wa u/ e osc . In the 
adiabatic regime, this ratio is small, and so 9 B <C 1- Then, as p precesses, the off-diagonal 
deviations are given by p X)V ~ |p| 6 B ~ (e sc)(e W aii/eosc) ~ ^waii- Meanwhile, the deviation in 
the ^-component giving the diagonal densities is 5p z ~ |p| (1 — cos9 B ) ~ e osc 9 B ~ e^/eosc- 
Thus we can neglect the feedback of 5p z to p XtV , and using the decoupling approximation to 
solve for j\ 2 = p x + ip y is justified. 

However, in the non-adiabatic regime, e wa n > e osc , the precession angle is order 1. Then 
the deviations in p x , y ,z are all the same order, 5p X)VtZ ~ |p| ~ e osc , and none of them can be 
neglected or approximated as being in equilibrium. The magnetic analogy makes clear why 
the decoupling approximation breaks down in the non-adiabatic regime — if the entire vector 
p precesses with a large angle, all components change simultaneously with equal magnitudes, 
and there is no notion of decoupled evolution of the individual components. 



C. Power Counting of Diagonal Solution 

For the diagonal densities, we can power count deviations of fa, fa from equilibrium as 



in Eq. (60). We find in the non-adiabatic regime, 5 fa ~ e^ all , so they certainly cannot be 
neglected. In the adiabatic regime, 8 fa ~ e^ all . Then it appears that the power counting of 
[9] is justified in this regime. However, this is not the case. 

First, deviations of diagonal distributions from equilibrium survive much farther into the 
unbroken phase (z > 0) than off-diagonals, as we found in Fig. |2j In Fig. [3] we found that 
the part surviving to very large z is well approximated by a chemical potential that survives 
even after kinetic equilibrium {fa = n B {ui — /ij), /i 2 = 0) is reached much earlier. This is 
because collisions drive to kinetic equilibrium on a length scale L co \ h but the chemical 
potential is damped away on a longer diffusion length scale L,as ~ L co \\/v w . Thus, for very 
large z, the diagonal 5 fa should not be neglected even in the adiabatic regime, although 
they are formally suppressed relative to /12 in the region of the wall. 



Second, as we argued in Sec. [TV] one is ultimately interested in the total integrated charge 
l2 B diffusing into the unbroken phase, and the amount of this charge is governed by the 
diffusion length Ldiff- in our model L^g is actually the largest length scale in the problem. 
Thus, although in the adiabatic regime the deviations of fa, fa from equilibrium are formally 
(9(e^ all ), their contributions to are order e^ all / (econf^,) due to the measure of integration, 
and thus parametrically larger than naively expected. In other words, even in the adiabatic 
regime, one should not neglect deviations from equilibrium in the diagonal densities. 



D. Numerical Comparisons and Diffusion Tail 



In Fig. |6| we compare our exact solution for /i2(k, z) and /i2(k, z) (solid lines) with the 
approximate solution Eq. (59 ) (dashed lines) in a typical adiabatic bin and non-adiabatic bin, 
given in terms of the components p x , y ,p X} y in the Bloch decomposition Eq. (45). Overall, the 
qualitative behavior is similar. Quantitatively, the solutions disagree in the non-adiabatic 
bins, while they are somewhat closer in the adiabatic bin, consistently with the power 
counting estimates given above. The remaining discrepancies may be due to effects missed 
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FIG. 6: Numerical results for the off-diagonal components p Xt y(k, cos z) of the particle density 
matrix (left panels) and anti-particle density matrix p x ,y(k, cos'dk, z) (right panels). Solid lines rep- 
resent solutions of the full equations. Dashed lines represent solutions of the approximate equations 
(51). The upper panels correspond to a typical adiabatic bin with k/T = 0.5, cos = 0.875, while 
the lower panels correspond to a typical non-adiabatic bin with k/T = 3.5, cos = 0.875. 



by the ansatz Eq. (56) for the collision term, on which we comment in Appendix Bl and due 



to neglecting the feedback of 5/11,22 to /12 through the wall-induced source in Eq. (53), on 
which we comment in Appendix [Cj 

The most significant discrepancy, however, between the treatment of Ref. [9] and ours 
arises in the diagonal distributions, and is illustrated in Fig. [7j Our exact solutions exhibit 
diffusion of the diagonal flavor densities n^^ deep into the unbroken phase (where they are 
equal to ^1,2), while this phenomenon is absent in the treatment of Ref. [9]. The reason is 
that the decoupling imposed in the source term in Eq. (55) does not allow the off-diagonal 
densities /i2,2i to source the diagonals /n,22- Thus, in Ref. [9J, the flavor-diagonal densities 
are obtained simply by rotating ^12,21 back to the flavor basis. But this generates nonzero 
ul^r only inside the wall, not outside, where the mass and flavor bases coincide. On the 
other hand we showed above that in a consistent power counting the fully-coupled source 
Eq. (53) must be kept. Once the diagonals fu 22 are sourced, through collisions they diffuse 



outside the wall into the unbroken phase. This is the main mechanism to generate in the 
unbroken phase. 
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FIG. 7: Charge density profiles ul(z) from the solution of the full equations (solid line) and the 
approximate decoupled equations ( dashed line ) that mimic the procedure of Ref. FJ|/. The left panel 
correspond to the off-resonance regime ttil/T = 2.6,mn/T = 2, while the right panel corresponds 
to the resonant regime rriL/T = 2.2, tur/T = 2. 



VI. CONCLUSIONS 

Electroweak baryogenesis is an attractive and testable explanation for the origin of the 
baryon asymmetry of the Universe. However, quantitative baryogenesis predictions for col- 
lider and low-energy precision measurements are obscured by orders-of-magnitude discrepan- 
cies between different treatments of the charge transport dynamics during the electroweak 
phase transition, due to different approximations. Since the final baryon asymmetry de- 
pends on how much CP-asymmetric charge is generated and transported into the phase of 
unbroken symmetry, a quantitative computation of these dynamics is essential. 

This work has provided an important step (beginning with Ref. [IB]) toward investigating 
these discrepancies and providing a quantitatively robust formalism for baryogenesis com- 
putations. Within a generalized gradient expansion, we derived Boltzmann equations for 
the particle and antiparticle density matrices for a two-flavor scalar system with an inhomo- 
geneous, CP-violating mass matrix. (One may regard the scalar species as a toy model for 
squarks in supersymmetric extensions of the Standard Model, with an inhomogeneous mass 
matrix arising from the spacetime- varying Higgs field during the phase transition.) Our 
Boltzmann equations are "flavored" in the sense that they keep track not only of the occu- 
pation numbers of individual states, but also of their coherence. These equations account for 
flavor oscillations in a non-homogeneous background and interactions with a thermal bath. 
In contrast to previous treatments, we have kept the full matrix structure of the collision 
term and have not resorted to the usual diffusion approximation. However, it is clear that 
diffusion does emerge from our full numerical solutions. From our analysis a very simple 
physical picture emerges: at leading order in gradients of the mass matrix, CP asymmetries 
arise from coherent flavor oscillations induced by spacetime-dependent mixing. 

By virtue of our simplified model, we solved the Boltzmann equations numerically without 
ansatz for the form of the density matrices. We illustrated several important physics points: 

• The largest departures from equilibrium and the largest contributions to the CP- 
violating asymmetries arise for states evolving non-adiabatically across the phase 
boundary, with momenta k satisfying L osc = 27rt> rcl (k)/(u;ik — w 2 k) > L w , where L w is 
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the bubble wall thickness (i.e. the length over which the off-diagonal elements of the 
mass matrix vary). 

• The CP-violating flavor-diagonal charge densities, generated by flavor oscillations 
within the wall, diffuse into the unbroken phase and are not localized near the wall. 

• The enhancement of charge generation in the non-adiabatic regime manifests itself as 
a resonance when w,l ~ tur. This "resonant regime" is governed by the condition 
V^l — mill ~ (i- e -' the width of the resonance is controlled by L w ). In MSSM- 
like models, L w ~ 20/T, so for mass differences as big as — ~ T/2 or about 
50 GeV, it becomes important to account consistently for modes that evolve non- 
adiabatically, as we have done. 

These findings demonstrate the crucial importance of keeping track of full coupled evolution 
of all components of the density matrix / to capture the dominant contributions to flavor- 
diagonal charge densities that diffuse into the unbroken phase. 

We also compared our results, within the context of our toy model, to the formalism 
of Ref. [9j. All previous baryogenesis computations have relied on an implicit or explicit 
decoupling of the dynamics of diagonal and off-diagonal densities, whereas our results here 
do not and achieve exact numerical solutions for the full density matrices. We have shown 
that the power counting of Ref. [H] (leading to decoupled equations for diagonal and off- 
diagonal densities) is inadequate in the non-adiabatic regime, in which the CP-violating 
effects are maximal. We have solved our simple model of mixing scalars according to the 
procedure outlined in Ref. [9], finding dramatic differences in the charge density profiles: 
the full solution displays significant diffusion into the unbroken phase that is absent in 
the approximate treatment. This difference can be directly traced to the approximation 
of decoupling the kinetic equations for diagonal and off-diagonal sources. This difference 
may have a potentially large impact on electroweak baryogenesis calculations. Our analysis 
indicates that the approach of Ref. [S] largely underestimates the CP-violating densities 
in the unbroken phase, which in turn might induce a large underestimation (by one order 
of magnitude or more) of the produced baryon asymmetry. Within the simple model of 
mixing scalars, however, we cannot address this in a quantitative manner, as we still need 
to introduce fermions. An additional CP- violating source, the semi-classical force [TTHT5] . 
can arise for fermions, but not scalars. The relative magnitude between this source and the 
resonant, mixing-induced source studied here remains an important open question. 

The resolution of current discrepancies and a more robust phenomenological analysis 
require the following additional steps, which are currently under investigation: 

• Transport equations for fermions with an inhomogeneous mass matrix along the same 
lines as our analysis of mixing scalars, including the resonant mixing-induced source, 
the semi-classical force, and elastic and inelastic scattering processes. 

• Identification of diffusion equations with appropriate oscillation-induced sources that 
correctly capture the physics of the full kinetic analysis. This should be a good de- 
scription of the system in the unbroken phase, where the mixing angle vanishes and 
flavor oscillations no longer occur. 

These developments, building upon the methods we have introduced here and in Ref. [T6] . 
will make possible rigorous and tractable predictions for charge transport in realistic sce- 
narios of EWB. 
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Appendix A: Collision terms 



We evaluate the interaction terms in the Boltzmann equations coming from Jz?i n t, given by 
Eq. (|2]). We work in the mass basis below, but omit the subscripts m for brevity. We assume 
that the A bosons (mass ttia) are in thermal equilibrium, with temperature T 3> m^. 

The self-energy functions II can be computed following a perturbative expansion in i/l,r, 
detailed in Ref. [16J. At linear order in yi,R, only the [II , G<] term receives a contribution, 
shown by Fig. (|8p,): 

a * (M ) = K (l) (£ + /|*_L). (Ai) 

The first term gives a thermal mass shift: m 2 L R — > m 2 L R + Ul,rT 2 /24. The second term 



(e p = y/p 2 + m A ) is the usual zero-temperature divergence that can be absorbed by renor- 
malization. 

At second order in Ul,r, scattering and annihilation processes arise from the imaginary 



part of Fig. (|8p) and are given by Eq. (32). The scattering term [(j)(k)A(p) -h- <f)(k')A(p')] is 
given by 

*tt /U =-^rJ T^Lo / f (2-) 4 S\k +p-k'-p') (A2) 



Ak° J (2n) 3 2k>° J (2n) 3 2e p J (27r) 3 2e p / 
x [{/(k), Y{1 + /(k'))y}/A(p)(l + /a(pO) - {(1 + /(k)), Yf(k')Y}(l + /a(p))/a( P / 

where /a(p) = ^s(ep) = l/( ex P( e p/^) — 1) is the distribution function the A bosons. The 
annihilation term [^(k^^k') <H- A(p)A(p')] is given by 

1 f d 3 k' f d 3 p f d 3 p' /0 _, 4r4 . 



W, fU =~ W J J J W W + V-P- V') (A3) 

[{/(k), Yf(k')Y }(1 + f A (p))(l + f A (p')) - {(1 + /(k)), F(l + fW)Y}f A (j,)f A <j,')] . 



X 



These expressions are strongly reminiscent of the corresponding single-flavor collision terms, 
except for their "non-abelian" structure: the distributions (i.e., density matrices) and scat- 
tering/annihilation matrix elements (Y) do not commute. The total collision term is 

W, /] = 9* (W, /Iscat + W, fUn ) , (A4) 

where we include the additional factor g* ~ 200 to mimic the true number of degrees of 
freedom in the plasma during the EWPT. 
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(a) (b) 

FIG. 8: Leading-order self-energy graphs that induce the collision terms in the Boltzmann equations, 
corresponding to (a) thermal mass corection from coherent foward scattering, and (b) non-forward 
scattering (4>A o <f)A) and annihilation ((fxjy <-> AA). 



There is one important subtlety: what is the dispersion relation for k° (and fc'°)? 16 
Working to linear order in e, it is consistent to set k° = ui^ in since corrections will be 
0(e co ne osc ). However, in that approximation, the collision term relaxes the density matrices 
to a "false equilibrium" /, / — » diag(ng(o;k), ris(aJk)), rather than the "true equilibrium" 
/,/—>■ diag(nB(wik), ^e(w2k)), with the correct dispersion relations. This distinction is 



crucial: according to Eq. (33), the CP- violating source term vanishes if /, / are proportional 
to the identity, which they are when "falsely" equilibriated. Hence, the C(e° sc ) collision 
term quenches charge generation, an unphysical effect. 

To remedy this issue, when evaluating the collision term we replace the lowest order 
solution to the constraint equations given in Eq. [24] with 



Gg(fc,z) = - m%) [e(k°)(5 tj + fij{k,x)) + ei-k^fiji-^x) } 

Gf J (k,x) = 2n5(k 2 -ml) [e{k )f ij {k,x) + 6{-k°){5 ij + f ij {-k,x)] . 

with mfj = l/2(mf + m?), implying 



(A5) 



k° = u& = \ Uk 1 ~ 3 . (A6) 
I Wk it 3 

This procedure can be viewed as the resummation of a class of 0(e co u X Cqsc) corrections to 
the collision terms, corresponding to dynamical effects over length scales L mfp /e™ sc (n > 1). 
It is physically justified to neglect these corrections for off-diagonal modes: these modes 
are damped to zero on a scale L m f p and are not sensitive to longer scales L m fp/e" sc . In 
contrast, diagonal modes do not equilibrate on a scale L m f p ; they approximately equilibrate 
to /, / ~ diag(ris(o>k), ^B(o>k)), but only reach true equilibrium over longer scales L m f p /e2 sc - 
Therefore, by adopting Eq. ( |A6[ ), we are evaluating collision terms involving diagonal modes 
to all orders in e osc , as is required to treat equilibration properly. Progress in evaluating ^ 
to all orders in e osc for both diagonal and off-diagonal modes has been made in Ref. j3(i| . 
With this prescription, a compact matrix expression for ^ is no longer possible and we 



16 This discussion applies to both particle (k° > 0) and antiparticle (fc° < 0) poles. In evaluating c <§ , we 
have k°, k'° > 0, where all negative energies are made positive through a change of variables. 
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must write all mass basis indices explicitly. The collision terms are 
f d 3 k' 

K™V, f] = 7^ (S ia Y bc Y dj + Y lc Y da 5 b3 ) (A7) 



,/ (2vr) 3 

x (^(k, k') (1 + /(k)) a6 f cd (k>) - i?S° Ut (k, k') / a6 (k) (1 + f(k')) cd 
C fi 3 k' 

KT n if, F\= Jj^ + Y lc Y da 5 bj ) (A8) 

x (« n (k, k') (1 + /(k)) a6 (1 + f(k')) cd - KZr\k, k') f ab (k) f cd (V) 
The scattering kernels are 

^r/ n (k, k') = TnB{t h o) . tf(t 2 - t 2 ) log ( 1 + (A9) 
abed v , ; Mirtugug V o; & V 1 + ^( t +)/ 

JST(k. k') = T(1 + ^ (to ] } £(t 2 - 1 2 ) log ( 1+n B<f-\ ) (A10) 
abed v , ; 647rto;f^ v oJ B \l + n B (t + )J v ; 

and the annihilation kernels are 

^ (k ' k } = 128^ « ^° " S " bg U(* + W(-* + ) J (AU) 

iCT (k, k') = T(1+ng fe (S ° ) d ) £( S 2 - , 2 - 4m 2 ) log ( n B^-)nB(-s.) \ 

abed V , ) USTTSutfug V ° ^ V™b(s+W-S+)/ 

where 

*=|k-k'|, to = a^-a^, t± = ±| + ^l + 4my(t2-^), (A13) 
s =|k + k'|, So = o;f + < s± = | ± | ^/l + 4my ( S 2 - a g) . (A14) 
From these expressions, one can verify several facts. First, detailed balance is satisfied since 

^r/ ut (k, ko = e ^ T i?r/ n (k, ko , icy ut (k, ko = e ^ T rtj u (k w . cms) 

Second, ^ vanishes for 

f(k) = ( nB ( Ulk ~ ^) ^ /(k) - /'^(^Ik + ^i) 



n B (w 2 k - A*2) / V ^B(w2k + /i2) 

(A16) 

with chemical potentials /ii^. (If V is diagonal in the mass basis, 01,2 charges are separately 
conserved; otherwise only total charge 4>i + 02 is conserved and = /i 2 .) Third, the 
continuity equation is satisfied provided Tr J <i 3 /c/(27r) 3 ( < ^'[/, /] — ^[/, /]) = 0, which follows 
from the relations 

^'° Ut (k, k') = i?r 6 ' in (k', k), iCT(k, k') = (k', k), (in out) . (A17) 
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Appendix B: Decoupling in the collision term 



In Ref. [Sj, the ansatz Eq. (56), = — T i2 f 12 , was made for the off-diagonal collision 



term. Here we consider corrections to this ansatz and when it may be justified. 



In general the collision terms given by Eqs. (A7) and (A8) have much more complicated 
structure than this simple ansatz. We can simplify them somewhat by working to linear 
order in deviations from equilibrium, taking 

fij = fij 1 + n B (uJij)(l + n B (uij))5fij , (Bl) 

where = n B (ui)5ij, and 5 fij ~ e for some small e. We have factored out n B {l + n B ) in the 
(9(e) term for later notational convenience (cf. [31 J ) . Then, we linearize the collision terms in 
5f. At O(e ), the collision terms vanish (as required for the equilibrium distributions). Now, 
consider the off-diagonal collision terms ^12 linearized in 5f. (We will study the structure 
of the diagonal collision terms in future work.) It can be organized into three sets of terms, 

C/P <v?loss 1 cfig&in . ^source /no"! 

75 12 — ©12 T ®12 ~r ®12 " 



The loss term takes the form of the ansatz Eq. (56) but with a k dependent relaxation rate, 
^(k^) = -T 12 (k,z)f 12 (k,z) = -T 12 (k,z)n B (ui 2 )(l + n B (u 12 ))5f 12 (k,z), (B3) 
where 

d 3 k' 



ri 2 (k,z) = —, r / 

1 + n B {u 12 ) J 



l + n B (u 12 ) J (2tt) 3 

(Y* + Y 12 Y 21 ) [R?£r\K k')(l + n B {u>l)) + R^iK k> fl (a&)] (B4) 



x 



+ (Y 2 \ + Y l2 Y 21 ) [RT 2 ^r\K k')(l + n B (u£)) + R^iK k')n B (u&)] 



Meanwhile, the "gain" term is 

d 3 k' 2 



^ 2 ain (k,,) = ^/^E 

x {R^\k^)n B {ul)[l + n B {uj 1 ^ 

- i?~ t (k,k)n B K)^(^)[r 11 r 22 5/ 1 2(k) + Y? 2 5f 21 (k')}} , 

and the "source" term is 

d 3 k' 



(B5) 



^r cc (Kz) = 9 *Y 12 J — J2 Y 3 
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x {i?^° ut (k,k')n B K)[l + n B (u3i,)][5f n (k') -Sf u (k)) 
- R-'° ut (k, k')n B (4)n B (4,)[54(k / ) + <5/«(k)] 



(B6) 
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so named since deviations of diagonal distributions of equilibrium act as a source for fa 
through a nonzero off-diagonal coupling Y 12 . 

The ansatz Eq. (56) misses the effects of both the collisional gain and source terms. In 
the limit of flavor-blind interactions, Vl = Vr = Vi we have Y 



12 



and Y 



11,22 



so the 



source terms vanish. However, a part of the gain term still remains. The remaining terms 
(so-called "noise terms" in [12]) are often neglected by assuming Sf(k!) to be a randomly 
fluctuating variable, causing the integral over k' to be suppressed relative to the loss term. 
Here a similar suppression may happen because of oscillations of 8 fa in both k' and z with 
frequency U\ — 0J2, but there is otherwise no a priori reason to drop these terms. 

In the comparison shown in Fig. [7] which uses the baseline parameters of Table [IJ it is the 
case that Y 12 < Y 112 2, but still nonzero. Thus the approximate solution using the ansatz 
Eq. (56) for the collision term misses the collision- induced source Eq. (B6) for fa. This 
is one likely cause of the smaller normalization of the approximate solution even in the 
adiabatic regime where Eq. (59) is otherwise valid. 



Appendix C: Power Counting the Off-Diagonal Distribution 



In Sec. VB we argued that the fa in Eq. (59) is actually 0(e° all ), not C?(e wa u) as argued 
in [9], and therefore could not be neglected in the source term Eq. (53). However, since in 



our power counting fa is still 0(e osc ) (see Eq. (58)), it is fair to ask why we still do not 
neglect the terms containing fa,2i m Eq. (53) since they become 0(e 2 ). The reason is that, 



by power counting this way, every component of the source term Eq. (53) becomes 0(e 2 ). 
The diagonal and off-diagonal components are all the same order and should all be kept. 
Counting consistently, the leading nontrivial terms in the kinetic equation become 0(e 2 ). 
Deviations of / from equilibrium can be counted as one power of e. Then the oscillation 
term is e osc e, the source term is e wa ne, and the collision term is e int e. So the derivative d z 
on the left-hand side of the kinetic equation always brings down at least one e when acting 
on /, and the whole kinetic equation begins (nontrivially) at 0(e 2 ). The upshot is that, 
counting consistently at this order, the full coupled matrix structure of the kinetic equation 
must be kept. 

This exercise in power counting also tells us that there is a regime in which the solution 



Eq. (59) from Ref. [9] is a relatively good approximation for fa- In the adiabatic regime, 
when L w 3> L osc , the factor exp(— iAuz") cuts off the range of integration to be of order 
L osc (over a larger region, oscillations average fa out to zero). Then fa in Eq. (59) is order 



fa ~ (range of integration) x £ 12 x [n B (ui) — n B {ui 2 ) 



x e wall x Cose — ^wall (CI) 



which is consistent with Ref. [9]. This is just the part of /i 2 sourced by the equilibrium 
diagonal distributions (call it fains])- One can show the deviations 5/11,22 of the diagonal 
distributions from equilibrium sourced by fains] is then order e wa u(e wa n/e osc ), which then 
feeds back to source an additional part of fa (call it fa[8f)) of order e wa n(e wa n/e osc ) 2 . Since 
e wa ii <C e osc in the adiabatic regime, these additional deviations are suppressed relative 
to fains]- So Eq. (59) should be a good approximation for fa in the adiabatic regime. 



However, as explained in Sec. |VB[ one should not decouple fa from the evolution of the 
diagonal densities in the source Eq. (53), even though 5/11,22 are formally 0(e^ all ) in this 



regime. Neglecting this source, and thereby neglecting deviations from equilibrium in fa, 
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precludes the existence of diffusion and therefore grossly underestimates the total charge in 
the unbroken phase. 
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